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the effect this degree of degeneracy has on the prospects for direct detection of relic nentralinos in 
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Introduction 

The perplexing question of the nature of the missing 
non-baryonic matter in the universe continues to be an 
active area of research. In supersymmetry, the lightest 
supersymmetric particle (LSP) tends to be a neutral gau- 
gino, which is stable if the assumption of R-parity conser¬ 
vation is imposed, and its weak-scale self-interaction rate 
in the early universe make it an attractive candidate to 
fill this role. Not surprisingly, then, relic nentralinos con¬ 
tinue to be the favored candidate for the missing matter 
density of the universe. 

Most phenomenological studies of the parameter space 
in supersymmetric models that give rise to the cosmo¬ 
logically preferred thermal relic density of 0.1 < 

< 0.3 focus on the minimal supergravity model 

(mSUGRA) with its unified gaugino mass mi /2 and uni¬ 
fied scalar mass mq. In such constrained models of the 
MSSM^ a typical point in the rno) plane that 

would otherwise be quite “natural” from the point-of- 
view of solving the hierarchy problem of the electroweak 
sector tends to produce far too much relic neutralino den¬ 
sity in the early universe to be compatible with the known 
age of the universe. A sizable portion of the parameter 
space that does not come into conflict with the bound 

< 0.3 is instead in conflict with the bound on the 
lightest CP-even Higgs mass coming from LEP, particu¬ 
larly for low values of the parameter tan/3. 

This disturbing situation was pointed out in our pre- 


^ The class of models which we refer to as mSUGRA is also com¬ 
monly known as the constrained MSSM, or “CMSSM,” in the 
literature. 


vious study [||, where it was observed that significant 
improvement can be achieved if the assumption of univer¬ 
sal gaugino masses is relaxed. The problematic behavior 
of the mSUGRA model is a result of the overwhelming 
B-ino content of the LSP and the large mass difference 
between the LSP and the next lightest supersymmetric 
particle (NLSP) which holds over most of the parame¬ 
ter space in such models. Combined with the need for 
a large universal gaugino mass to accommodate a Higgs 
mass of 115 GeV these properties lead to an excess of 
relic B-inos outside of certain special areas where reso¬ 
nant annihilation or coannihilation processes dominate. 
When this assumption of gaugino mass universality is 
not valid - in particular when the W-ino mass M 2 is less 
than the B-ino mass Mi at the initial supersymmetry 
breaking scale Ayy - the increased W-ino content of the 
LSP and reduced mass difference Am = m^± — m^o be- 
tween the LSP neutralino Xi and NLSP chargino can 
have dramatic effects on the depletion of relic neutrali- 
nos in the early universe. This will open up new areas of 
MSSM parameter space with acceptable relic densities. 
Quite apart from this improvement in the dark matter 
arena, such nonuniversalities in the gaugino sector are 
a common feature of string-derived effective theories. It 
was further suggested recently j|] that nonuniversalities 
in the gaugino sector may also play an important role 
in ameliorating the apparent fine-tuning of the MSSM 
required in the electroweak sector as well. 

In the current work we utilize newly-available calcula¬ 
tion tools to extend the analysis to include the important 
large-tan/3 region which has been of great interest re¬ 
cently in the context of minimal supergravity models. We 
pay particular attention to the behavior of that area of 
parameter space in which resonant annihilation through 
heavy Higgs states occurs when gaugino mass nonuniver- 
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sality is introduced. We find that the viable parameter 
space is enlarged, with significant regions allowing for 
heavy scalars and relatively light gauginos. The physical 
significance behind this new parameter region is investi¬ 
gated and the sensitivity of these results to the value of 
the gluino mass is explored. This analysis forms the bulk 
of Section 

Just as changing the relative size of the gaugino mass 
parameters can have large effects on the number and na¬ 
ture of coannihilation channels available to the relic LSP 
in the early universe, so too can it have profound im¬ 
pacts on the interaction rates of surviving relics with cur¬ 
rent dark matter detectors. In Section || we investigate 
the prospects for direction detection of relic neutralinos 
in our local halo in the next round of future detectors. 
While nonuniversalities that give rise to a relic density in 
the preferred range typically generate interaction cross- 
sections too small for current detectors such as DAMA 
or CDMS, we find a significant region of the preferred 
parameter space will be probed by the next generation 
of Germanium and Xenon-based detectors. 


Lastly, in an effort to quantify the difference between 
the mSUGRA paradigm and the case with nonuniversal 
gaugino masses, we adopt a fine-tuning sensitivity pa¬ 
rameter similar to that of Ellis and Olive Q . We calcu¬ 
late the sensitivity of the neutralino relic density to small 
variations in the high scale parameters, including input 
values such as the top and bottom quark masses. In 
Section III we compare the typical amounts of tuning in 
universal and nonuniversal scenarios, and comment upon 
the implications, before concluding. 


I. RELIC DENSITY OF LSP NEUTRALINOS 


We are interested in studying the present density of 
neutralinos assuming they were once in thermal equilib¬ 
rium with radiation in the early universe. The lightest 
such neutralino will be the LSP in all models we will con¬ 
sider. The mass and couplings of the neutralino sector 
at tree level are determined by the mass matrix 


/Ml 0 

0 Ma 

— sin 01 V cos/JM^ cos 9w cos PMz 
\ sin dw sin PMz — cos 9w sin PMz 


— sin 0 iv cos PMz 
cos 9w cos PMz 


0 


sin 9w sin /3Mz 
— cos 9w sin PMz 


0 


( 1 ) 


which is written in the {B,W, H^, H^) basis, where B 
represents the B-ino, W represents the neutral W-ino 
and and are the down-type and up-type Higgsi- 
nos, respectively. The parameters Mi and Ma in are 
the soft supersymmetry-breaking gaugino masses for the 
B-ino and W-ino, respectively while /r is the supersym¬ 
metric Higgsino mass parameter. The ratio of Higgs vevs 
is given by the parameter tan/3 = Vujvd with the vac¬ 
uum expectation value of the up-type Higgs and Vd the 
vacuum expectation value of the down-type Higgs. We 
will assume that all of the parameters in ( ^ are real for 
this analysis.^ 

The content of the lightest neutralino can be parame¬ 
terized by writing the LSP as 


X° = /ViiB + iVi2W + Ai3MS + iVi4M°, (2) 

with the normalization iVfi-I-iVf 2 -I-A^i 3 -I-iVi 4 = 1. When 
the parameter /i is much larger than the Z-boson mass 
it is clear that the proportion of W-ino content to B- 
ino content in the LSP is determined by the ratio of 
their soft masses, which we denote by r = M 2 /Mi. In 


^ See, however, Ref. ^ for the implications on neutralino relic den¬ 
sities when this assumption is relaxed. 


mSUGRA this parameter is defined to be r = 1 at some 
high-energy input scale, usually taken to be the GUT 
scale Auv = A^u^ = 2 x 10^® GeV. Upon renormalization 
group (RG) evolution of these parameters to the elec- 
troweak scale this implies Mi ~ 5 M 2 and the lightest 
eigenvalue of (|^ is thus overwhelmingly B-ino in content 
{Nil — 1). In a more general scenario the neutralino 
sector would be defined by four parameters at the elec- 
troweak scale: {tan/3, /r, r. Mi}. 

In practice the magnitude of the /it parameter is deter¬ 
mined at the electroweak scale by enforcing proper elec- 
troweak symmetry breaking (EWSB). Here this will be 
performed by computing the complete one loop corrected 
effective potential Vi_ioop = Utree + AUad- The correc¬ 
tion AUrad includes radiative effects from the entire su¬ 
perpartner spectrum and thus involves all of the param¬ 
eters of the soft supersymmetry-breaking Lagrangian. 
The effective /r-term p, is calculated from the EWSB con¬ 
dition 


r,2 _ 

H- — 


- {r 


"'Hu 


5rn\j^ 


tarn /3 — I 


) tan /3 1 „ 

^- 2 ^- 


(3) 


where Stuh^ and Sttih^ are the second derivatives of the 
radiative corrections AUad with respect to the up-type 
and down-type Higgs scalar fields, respectively. The sign 
of the ^ parameter is not determined and left as a free 
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FIG. 1: Neutralino Relic Densities for tan/3 = 5. The 
preferred neutralino relic density regions (light shading) are 
displayed for r = 1 (left panel) and for r — 0.6 (right panel). 
Dark shaded regions in the lower right are excluded due to a 
charged LSP. Medium shaded regions on the left are excluded 
by LEP bounds on the chargino mass. The Higgs mass con¬ 
tour of irth = 113 GeV is given by the dotted vertical line. 
The dashed contours are a rough estimate of the direct detec¬ 
tion reach of the GENIUS (“GEN”) and XENON (“XEN”) 
experiments. 


FIG. 2: Neutralino Relic Densities for tan j3 = 35. The 
preferred neutralino relic density regions (light shading) are 
displayed for r = 1 (left panel) and for r — 0.6 (right panel). 
Dark shaded regions in the lower right are excluded due to a 
charged LSP. Medium shaded regions on the left are excluded 
by LEP bounds on the chargino mass. The Higgs mass con¬ 
tour of ruh = 113 GeV is given by the dotted vertical line. 
The dashed contours are a rough estimate of the direct detec¬ 
tion reach of the GENIUS (“GEN”) and XENON (“XEN”) 
experiments. 


parameter. 

Determination of /r through (||), as well as the deter¬ 
mination of the one loop radiative corrections to the neu¬ 
tralino mass matrix (Q), involve the remainder of the su¬ 
persymmetry breaking soft terms. These soft terms will 
also be necessary to translate low energy values in (|^) 
into high energy input values by solving the renormal¬ 
ization group equations (RGBs). We wish to isolate the 
effects of gaugino mass nonuniversality so we choose to 
consider a universal scalar mass toq and universal trilin- 
ear A-term Aq as in minimal supergravity. We are thus 
led to consider the seven-dimensional parameter space we 
refer to as rSUGRA given by 

{tan/3, sgn(/i), r, Mi, M 3 , mo, Aq} . (4) 

The effect of the trilinear coupling Aq is negligible for 
the physics we are considering here so we will set Aq = 0 
throughout. We will also restrict our attention to positive 
/i. This is very similar to the parameter space originally 
studied by Mizuta, Ng and Yamaguchi |^, where a trio 
of specific values for the parameter r were investigated. 
The effect of gaugino mass nonuniversality on neutralino 
dark matter in certain specific models has also been com¬ 
puted recently in Ref. || and The latter appeared as we 
completed this analysis and we will comment on areas of 
agreement with that work during the following sections. 

Translation of the parameter set (Q) from their values 
at the high scale (we will assume Auv = Aqut) is per¬ 
formed using the Fortran code SuSpect which solves 
the RGBs at two loops between the scale A^ux and Mz 
through an iterative procedure that takes into account 
proper thresholding of superpartners. Proper BWSB is 
imposed and the one loop corrected /r parameter is de¬ 
termined. The complete superpartner spectrum is calcu¬ 
lated including one loop radiative corrections. 


The soft supersymmetry breaking parameters at the 
weak scale are then passed to the G code micrOMEGAs |l^ 
to perform the relic density calculation. This consists of 
computing a thermally averaged annihilation cross sec¬ 
tion times LSP velocity and freeze-out tempera¬ 

ture xf = Tp/rn^o^ for the LSP Xi @1 l^li |^j @- 
Freeze-out occurs when the relic density of the LSP no 
longer tracks the density of relativistic degrees of freedom 
in the universe but is instead nearly constant. The freeze- 
out temperature is defined in micrOMEGAs as the temper¬ 
ature at which the relic density of Xi becomes roughly 
twice the value that would obtain in equilibrium. The 
power of micrOMEGAs lies in its inclusion not just of di¬ 
rect annihilation processes involving two neutralinos such 
as XiXi ^ but also of all possible coannihilation 

channels involving the relic LSP such as XiX^ ^ W^Z 
and Xi'? JT- These processes are known to be crucial 
to a proper calculation of the relic density whenever the 
masses of the LSP and of the coannihila ting pa rticle are 
nearly degenerate 


A. The case of gluino/b-ino equality 

As mentioned in the introduction, the typical 
mSUGRA analysis has the luxury of displaying results 
in the {mi/ 2 , mo} plane for a given value of tan/3. Our 
parameter set given in (^) introduces two new degrees 
of freedom for the gaugino sector, so we will begin our 
analysis by exhibiting the neutralino relic density in the 
(Ml, Too) plane for specific values of r and tan/3. In 
order to make contact with the familiar results of recent 
mSUGRA studies |2|, |^, |2|] we will set M 3 = Mi 

at the initial GUT scale in this subsection. Since the bulk 
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FIG. 3: Neutralino Relic Densities for tan j3 = 50. The 
preferred neutralino relic density regions (light shading) are 
displayed for r = 1 (left panel) and for r = 0.6 (right panel). 
Dark shaded regions in the lower right are excluded due to 
a charged LSP, while those in the upper left are excluded 
due to improper EWSB. Medium shaded regions in the lower 
left are excluded by LEP bounds on the chargino mass. The 
Higgs mass contour of mh = 113 GeV is given by the dotted 
vertical line. The dashed contours are a rough estimate of the 
direct detection reach of the GENIUS (“GEN”) and XENON 
(“XEN”) experiments. 


of the RGE effects of a given value of 1111/2 in unified mod¬ 
els is related to the value of M 3 , this choice most closely 
mimics the RGE behavior of the universal mSUGRA sce¬ 
nario. 

The results of this analysis are shown in Figures 01 
and H for values of tan/3 of 5, 35 and 50 respectively. 
In each of these figures, we display available parameter 
space in the (Mi, mo) plane, with Mi and mg being the 
soft supersymmetry breaking parameters at the initial 
(GUT) scale, for r = 1 and r = 0.6. The former value 
recovers the case of minimal supergravity, while the lat¬ 
ter value was observed in Ref. ^ to be the area where 
nonuniversality in the gaugino mass sector reaches a crit¬ 
ical value for which the preferred relic density becomes 
roughly independent of scalar mass values. The lightly 
shaded area in each region is the area of the parameter 
space preferred by cosmology, with 0.1 < < 0.3. 

The left-most contour (green) contour bordering this 
region is the contour of ~ ^.l while the right¬ 

most (red) contour bordering the region is the contour 
of = 0.3. Recent analysis of all applicable cos¬ 

mological data |3^ suggests a narrower region for lower 
values of fl^h^ may be preferred by observational data. 
We therefore provide the (blue) contour for fl^h^ = 0.2, 
within the preferred region, as a guide to this more re¬ 
strictive parameter space. 

The darkly shaded areas in each plot are ruled out by 
requiring the lightest neutralino to be the lightest super¬ 
partner (the region in the lower right) or by requiring that 
correct electroweak symmetry breaking occurs (the re¬ 
gion in the upper left on some plots). While it has become 
fashionable for dark matter calculational tools to include 
computations of indirect constraints on supersymmetry 
such as BR (6 —> sy) and 5(a^)|susY, we have chosen not to 


include such information in our figures to reduce clutter.^ 
We have checked, however, that with our choice of /r > 0 
the conservative l-cr lower limit on the branching ratio for 
inclusive decays BR (6 —> sy) > 2.2 x 10“"'^ coming from 
CLEG is satisfied for all of the parameter choices 
we exhibit. So too is the conservative upper bound |3|] 
on the SUSY contribution to the muon anomalous mag¬ 
netic moment (5(a^)|susY < 90 x 10“^° [Q for all but the 
slimmest region in the lower left corner of these figures 
- a region irrelevant to the parameter region we are ex¬ 
ploring in this paper. Two key experimental constraints 
which we do display are the Higgs and chargino mass lim¬ 
its from LEP. Regions excluded by the LEP limit on the 
mass of the lightest chargino to^± > 103 GeV |^, |3|] 
are indicated by the medium (magenta) shaded regions 
on the left of the plots. We have chosen to be conser¬ 
vative and indicate the contour of rrih = 113 GeV for 
the lightest GP-even Higgs mass by the dotted line as an 
indication of the area excluded by LEP 

Finally, the estimated detection reach of the proposed 
Germanium detector GENIUS, and the next-generation 
liquid Xenon-based XENON detector, are indicated by 
the heavy dashed lines labeled “GEN” and “XEN,” re¬ 
spectively. These are obtained by taking the estimated 
minimal cross section for neutralino-nucleon elastic scat¬ 
tering, cr^p, capable of producing a signal in GENIUS [stI] 
and XENON |Q, as a function of the LSP mass m^o. 
This is translated into a contour in the (Mi, mo) plane 
and normalized for a neutralino relic density of H^h^ = 
0.1. We will discuss these contours in more detail in Sec¬ 
tion 1 ^. 

We begin with the cases for low to moderate tan/3 in 
Figures ^and ^ which display similar features. For r = 1 
there is very little parameter space available for which the 
Higgs exclusion bound is not violated, particularly for low 
tan /3. One is required to be in the thin coannihilation tail 
where the lightest stau is almost exactly degenerate with 
the lightest neutralino. This situation is well-known in 
mSUGRA (the r —> 1 limit of rSUGRA). However, once 
r is lowered towards the critical value of r = 0.6 we see 
that a wide plume opens up, independent of toq, that 
is perfectly consistent with all experimental constraints 
arising from accelerators. For both low and moderate 
tan (3 this preferred region for r = 0.6 arises at relatively 
large values of Mi , implying relatively large values of the 
gluino mass, and thus occurs in areas of the parameter 
space where the Higgs mass can easily exceed the LEP 
limit. 

For extremely large values of tan/3, as Figure |3| demon¬ 
strates, the Higgs mass constraint is satisfied throughout 
most of the parameter space. Here resonant annihila¬ 
tion of neutralinos can proceed through an s-channel ex¬ 
change of neutral heavy Higgs states, particularly the 


® It should be noted that micrOMEGAs also computes such quanti¬ 
ties, as well as applying direct collider limits. 
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pseudoscalar. This generates an increase in the cosmo- 
logically preferred region, even for r = 1. The large size 
of this region is a direct result of the large width of the 
CP-odd Higgs mass in this high tan/3 regime and has 
generated much recent interest @ ® E3. il- 

The position, and indeed even the existence, of this 
large region for any given value of tan (3 is very sen¬ 
sitive to the Yukawa couplings of the bottom and top 
quarks |^. The pole masses we used in SuSpect for this 
analysis were = 175.0 GeV and = 4.9 GeV. 
This choice is slightly higher than the default setting of 
= 4.6 GeV in SuSpect, but we find this value al¬ 
lows us to match our results most closely with those of 
Refs. 


40, O for tan/3 = 50. This should be seen as 


evidence of the sensitivity of high tan /3 results to the sub¬ 
tleties of one’s RGE code and how one treats the NLO 
corrections to the Higgs sector."^ From the right panel 
of Figure ^ it is clear that the presence of this Higgs 
pole region continues to distort the allowed parameter 
space as we increase the degree of nonuniversality in the 
B-ino/W-ino system, but the sensitivity of this allowed 
space to Yukawa couplings is greatly reduced. We will 
return to this issue in Section HI below. 


It would be tempting to conclude from Figures 
through 1^ that relaxing the unification constraint on the 
gaugino masses will increase the supersymmetric param¬ 
eter space that can properly account for the cold dark 
matter of the universe. And while it is certainly true 
that any particular point in the (Mi, toq) plane - par¬ 
ticularly those at larger values of these soft parameters 
- can now always accommodate the cosmologically pre¬ 
ferred relic density for some particular choices of r and 
tan /3, the size of the allowed areas in these figures is 
somewhat deceptive. More information can be obtained 
by exchanging the variable mg for r, since it is this pa¬ 
rameter that is more directly tied to the physics that 
determines the eventual LSP relic density. 


For this reason we plot the neutralino relic density in 
the (Ml, r) plane in Figure ^ for tan/3 = 5 (the top pair 
of plots), tan/3 = 35 (the middle pair of plots) and for 
tan/3 = 50 (the bottom pair of plots). In this case we 
have fixed the scalar mass to two different values to iso¬ 
late the effects of the gaugino sector parameters alone. 
Displayed in this manner it is clear that the allowed re¬ 
gion in gaugino mass space is a narrow band that extends 
from the mSUGRA limit of r = 1 in the top left, where 
the region of direct annihilation between very light LSPs 
is ruled out by the LFP chargino mass limits, eventu¬ 
ally extending to lower values of r, and then connecting 
continuously with the stau coannihilation region for low 
values of toq. We believe that this space is the natural 
one for discussing neutralino dark matter in nonuniver- 
sal gaugino mass models: it conveys far more information 


For an analysis of mSUGRA parameter spar e i nvolving the de¬ 
fault SuSpect settings, one may consult Ref. |28[ 
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FIG. 4: Preferred Region in Gaugino Mass Space. Top 

panels are for tan/3 = 5, middle panels are for tan/3 = 35 
and bottom panels are for tan/3 = 50. The preferred neu¬ 
tralino relic density regions (light shading) are displayed for 
universal scalar mass mo = 500 GeV (left panels) and for 
mo = 1500 GeV (right panels). The dark shaded region in the 
upper right of the left panels is excluded due to the charged 
LSP, while the region on the left in the tan /3 = 50 case is ex¬ 
cluded by improper EWSB. The medium shaded regions on 
the left are excluded by LEP bounds on the chargino mass. 
The Higgs mass contour of nih = 113 GeV is given by the dot¬ 
ted vertical line. Contours of constant degeneracy parameter 
A are given by the thin near-horizontal lines. We have indi¬ 
cated by the dashed contours a rough estimate of the direct 
detection reach of GENIUS (“GEN”) and XENON (“XEN”), 
where applicable. 


on neutralino dark matter than the typical ™o) 

plane in that the correlation between degeneracy and 
relic density is readily apparent. 

To elaborate, we have plotted contours of constant de- 
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FIG. 5: Preferred Mass Degeneracy in the Gaugino 
System when M3 = Mi and tan/3 = 35. The left vertical 
axis measures the mass (solid lines) of the lightest and next- 
to-lightest neutralinos Xi and X 2 , respectively, as a function of 
the parameter r. The right vertical axis measures the degree 
of degeneracy (dashed line) between the LSP and the chargino 
as given by the parameter A, in percent, as a function of 
the same parameter. The cosmologically interesting region 
0.65 ^ r ^ 0.60 is indicated by the red dotted lines. 


neutralino are plotted as a function of the ratio r when 
Ml is held fixed at Mi = 500 GeV and toq is set to 
Too = 1000 GeV. The level crossing between a predomi¬ 
nantly B-ino LSP to a predominantly W-ino LSP occurs 
abruptly at r ci 0.55. But just prior to this transition, 
near r ~ 0.65 the degeneracy parameter A between the 
chargino and the LSP drops rapidly through the cosmo¬ 
logically preferred region, as is measured by the right 
vertical axis. Hence the justification, a posteriori, for 
focusing on the case r = 0.6 in Figures through 

Of course, the nature of the gaugino mass system is 
also determined by the value of the /r parameter at the 
electroweak scale - and in our procedure of fixing that 
parameter via the EWSB conditions this implies some 
dependence on the gluino mass M 3 and the various sca¬ 
lar masses through RG running. We emphasize that the 
importance of the gluino mass on the relic densities is 
of a secondary nature, through its indirect RGE effects. 
Nevertheless, these effects can be significant in certain 
regions of parameter space. 


generacy degree A defined by 


A = 


m+ — m^o 

Xi _^ 

rfi 


(5) 


where m is the average of the chargino and LSP masses. 
For moderate to large values of the LSP mass, as de¬ 
termined by the parameter Mi at the GUT scale, the 
preferred region of parameter space is in the neighbor¬ 
hood of r ~ 0.6 where the degree of degeneracy between 
the two mass eigenstates is roughly A = 5%. In this re¬ 
gion of the parameter space the degree of degeneracy is 
just sufficient to allow coannihilation processes between 
the chargino and LSP to sufficiently reduce the otherwise 
large thermal relic density of the lightest neutralino [||. 
Clearly, as this degeneracy is increased still further at 
lower values of the parameter r the coannihilation pro¬ 
cesses will be too efficient and the lightest neutralino 
ceases to be an attractive dark matter candidate with¬ 
out some method of non-thermal production. This is, 
for example, the case in models of anomaly-mediation 
for which Tamsb = M 2 /M 1 = 5/33 ~ 0.15 at the GUT 
scale Q. 

The importance of this special value of A, and hence of 
the phenomenological parameter r = M 2 /Mi, is no g reat 
mystery. As was pointed out by Griest & Seckel the 
coannihilating particle must be thermally accessible to 
the LSP near the time of freeze-out for coannihilation 
processes to have a significant impact on the overall relic 
density. Typical freeze-out temperatures in the parame¬ 
ter space we consider are on the order of Tp ~ m^o/20, 
implying a degeneracy on the order of a few percent is 
required to make coannihilation processes relevant. 

That the relic density should change so rapidly as r 
is diminished for fixed Mi is also not surprising when 
one considers the neutralino and chargino mass system. 
In Figure || the mass of the lightest and next-to-lightest 


B. The case of gluino/w-ino equality 


In Section ^ we looked exclusively at the case where 
M 3 = Ml at the GUT scale in order to reproduce the 
typical RGE behavior of universal models. In mSUGRA, 
with its unified gaugino masses, most of the physics we 
associate with a given value of TOi /2 is really the result 
of the RGE effects of the gluino. When studying neu¬ 
tralino relic densities the value of Mi at the GUT scale 
is important in its own right, since the LSP mass is typ¬ 
ically m^o ~ 0.4Mi in such scenarios. The separate im¬ 
pacts of tliese two different mass scales are obscured when 
Ml = M 3 , though this allows us to compare results to 
that of mSUGRA more readily. 


Rather than treat the gluino mass as a completely free 
variable we prefer to keep the analysis tractable by in¬ 
voking the relation M 3 = M 2 in this subsection. Such 
a scenario will lead to gluinos which are far lighter than 
in the previous section, especially in the area around the 
critical region of r ~ 0.6. What’s more, we believe such 
a relation is likely to be closer to the truth if gaugino 
mass nonuniversality arises in models based on heterotic 
string theory. 


For example, should gaugino mass nonuniversality 
arise at tree level in weakly-coupled string models it 
would probably do so via nonuniversal affine levels ka 
for the Standard Model gauge groups Qa, since gau¬ 
gino masses at tree level are determined by the auxil¬ 
iary field of the dilaton which couples universally to the 
gauge groups. Standard constructions for string mod¬ 
els enforce k = 1 for non-abelian gauge groups such as 
SU{3) and SU{2), but generally have no such restriction 
on the affine factor ky associated with the U{1) of hy- 
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FIG. 6 : Preferred Mass Degeneracy in the Gaugino 
System when M3 = M2 and tan/3 = 35. The left vertical 
axis measures the mass (solid lines) of the lightest and next- 
to-lightest neutralinos Xi and X 2 j respectively, as a function of 
the parameter r. The right vertical axis measures the degree 
of degeneracy (dashed line) between the LSP and the chargino 
as given by the parameter A, in percent, as a function of 
the same parameter. The cosmologically interesting region 
0.65 ^ r ^ 0.60 is indicated by the red dotted lines. 
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FIG. 7: Neutralino Relic Densities for tan/3 = 5 and 
M 3 = M 2 . The preferred neutralino relic density region (light 
shading) is displayed for r = 0.6. The shaded regions and 
contours are the same as those of Figure 0. 
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percharge.® Universality at tree level is recovered when 
the phenomenological choice ky = 5/3, corresponding to 
the standard GUT normalization of the hypercharge in 
the Standard Model, is made. But this assumption is 
typically not borne out in explicit string compactihca- 
tions where exotic hypercharge normalizations are quite 

. While such nonstandard hy- 
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common 

percharge normalizations make gauge coupling unifica¬ 
tion less straightforward, the presence of additional states 
in the theory with exotic hypercharges can reconcile the 
Standard Model gauge couplings with the universal cou¬ 
pling of the underlying string theory . 


A perhaps more palatable alternative for nonuniver¬ 
salities is that the tree level gaugino masses are indeed 
universal, but somehow suppressed through the mecha¬ 
nism of supersymmetry breaking and transmission to the 
observable sector. Such a situation is, in fact, quite com¬ 
mon in realistic string theories and can arise when the 
dilaton does not participate in supersymmetry breaking 
at all, or when its lowest component is stabilized using 
stringy nonperturbative contributions to its Kahler po¬ 
tential In both cases the resulting gaugino masses 
are loop-suppressed, implying that non-universal correc¬ 
tions are relevant. One class of terms arising at the one- 
loop level are those associated with the superconformal 
anomaly SH which tend to increase the gaugino mass 
Ml far more than that of M 2 and M 3 . 


In such a regime where M 2 « M 3 the level crossing that 
occurs in the neutralino system occurs far less abruptly. 
In Figure ^ we again show the masses of the LSP and 
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See, for example, the discussion of this point in Ref. and the 
references contained therein. 


NLSP neutralinos for Mi = 500 GeV, toq = 1000 GeV 
and tan /3 = 35. The gluino mass at the GUT scale ranges 
from 500 GeV at r = 1 to 125 GeV at r = 0.25 on the 
right, leading to much smaller values of the resulting /r 
term than in Figure and a smoother transition from 
B-ino like to W-ino like LSP. Nonetheless, the critical 
region where A ~ 5% continues to fall in the area near 
r = 0 . 6 . 

The effect of the lighter gluino mass is seen most dra¬ 
matically in the plots analogous to those of Figure |^. 
We return to the (Mi, toq) plane for tan/3 = 5 in Fig¬ 
ure^, where now we have the relation M 3 = M 2 = O. 6 M 1 . 
As in the right panel of Figure the preferred cosmo¬ 
logical region is centered around Mi ~ 1 TeV, though 
this region is narrower than in the previous case. More 
significant is the region of rapid neutralino annihilation, 
and neutralino-chargino coannihilation, through heavy 
Higgs states that now appears even at low values of tan /3. 
These heavy Higgs eigenstates are lighter than in a com¬ 
parable unihed model because of the diminished effect of 
the gluino mass on the ultimate size of the /i parameter, 
resulting in a smaller value than in mSUGRA. In the case 
of Figure ^ two different resonant regions can be resolved: 
the lower pole region being chargino-neutralino resonant 
coannihilation through charged Higgs states, while the 
upper pole region is neutralino-neutralino resonant coan¬ 
nihilation through the heavy neutral Higgs states H and 
A. 

As mentioned earlier, it is well known that the size 
and location in the {rni/ 2 , tan/3) space of this rapid 
annihilation region is a strong function of the top and 
bottom quark pole masses one assumes, as well as the 
treatment of the relevant corrections to the heavy Higgs 
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FIG. 8 : Neutralino Relic Densities for tan/3 = 35, 50 
and M 3 = M 2 ■ The preferred neutralino relic density regions 
(light shading) are displayed for r — 0.6 and tan /3 = 35 (left 
panel) and tan /3 = 50 (right panel). The dark shaded regions 
in the lower right are excluded due to a charged LSP, while 
those in the upper left are excluded due to improper EWSB. 
The medium shaded regions on the left-most edges are ex¬ 
cluded by LEP bounds on the chargino mass. The Higgs 
mass contour of rrih = 113 GeV is given by the dotted ver¬ 
tical line. The dashed contours are a rough estimate of the 
direct detection reach of the GENIUS (“GEN”) and XENON 
(“XEN”) experiments. 


masses.® The analysis of this section suggests that on 
top of this sensitivity one should add the sensitivity to 
the gluino mass, whose effects are felt through the RGBs. 

This effect becomes even more pronounced at higher 
values of tan/3, as in Figure where the heavier Higgs 
boson states become lighter still. For the case of tan/3 = 
35 in Figure ^ (the left panel) the resonant annihilation 
region has become quite pronounced as the region where 
2mj^o Ri iTiA, ttih moves to lower values of Mi. The size 
of the depletion region has also increased substantially 
over the case where M 3 = Mi and tan /3 = 50. This is the 
result of resonant annihilation through the pseudoscalar 
Higgs being possible at lower values of tan /3, where the 
width is narrower and the resonance annihilation ef¬ 
fects are more pronounced. In the center of this resonant 
annihilation region the neutralino relic density is very 
low, implying a much lower count rate for GBNIUS and 
XBNON and hence a lower direct detection reach in this 
area. By the time tan /3 = 50 is reached (the right panel 
in Figure |^) this resonant annihilation region no longer 
impacts the preferred relic density region and is confined 
(for this value of r = 0.6) to the region Mi ^ 1 TeV. The 
ultimate conclusion one should draw is that the presence 
of gluinos lighter than those in the universal case will gen- 


In fact, even the choices of as(Mz) and the exact definition of 
the GUT scale Auv = Agut can have effects on the pseudoscalar 
mass large enough to effect the position of this depletion region. 
We have chosen the value as{Mz) = 0.119 and determine the 
high scale Auv by the scale at which gi and 92 unify, usually 
close to the value of Agut = 1.9 X 10^® GeV. 
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FIG. 9: Preferred Region in Gaugino Mass Space for 

Ms = M 2 . Top panels are for tan/3 = 5, middle panels are 
for tan/3 = 35 and bottom panels are for tan/3 = 50. The 
preferred neutralino relic density regions (light shading) are 
displayed for nniversal scalar mass mo = 500 GeV (left pan¬ 
els) and for mo = 1500 GeV (right panels). The contours in 
this figure are identical to those of Figure W 


erally imply heavy Higgs resonant annihilation at lower 
values of tan/3 than in mSUGRA, as well as shifting this 
region to higher values of the universal scalar mass for 
fixed Ml and tan /3. 

It should be pointed out that by choosing to look in 
the (Ml, Too) plane for the case where r = 0.6 we have 
found that the preferred region in that plane tends to 
fall around Mi = 0{1 TeV) for both choices of M 3 
and all three values of tan/3. But this is an artifact 
of the choice r = 0.6. The LSP mass can be arbi¬ 
trarily light, up to the constraint imposed by the LFP 
bound on the chargino, by considering other values of 
r for which chargino-neutralino coannihilations are less 
efficient. This can be seen in Figure P, the analog to 
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Figure shown previously, where lighter LSP neutrali- 
nos can be obtained provided 0.7 ^ r ^ 0.9, depending 
on the value of tan /3. In this region the chargino and 
neutralino are farther apart in mass so that A ^ 20% 
and coannihilation processes are losing their relevance 
in favor of the standard annihilation processes of typi¬ 
cal mSUGRA scenarios. This region of light charginos 
is also the region most accessible to GENIUS and might 
be favored on general fine-tuning arguments related to 
electroweak symmetry breaking 

The distortions in Figure^, evident for low mo in the 
tan P = 5 case and for both mg values in the higher tan /3 
cases, are a result of the resonant annihilation region. In 
this region the correlation between chargino/neutralino 
degeneracy, as measured by A, and the neutralino relic 
density is destroyed. Note also that the reduction in relic 
density makes this region generally inaccessible experi¬ 
mentally, even by the proposed XENON detector. When 
this distortion is factored out the detectability reach in 
Ml of relic neutralinos in our local halo improves in the 
M 3 = M 2 regime over the case where M 3 = Mi. 


II. PROSPECTS FOR DIRECT DETECTION OF 
RELIC NEUTRALINOS 


We now wish to discuss the issue of direct detection 
of LSP neutralinos in more detail, couching the discus¬ 
sion in terms of the parameter space reach of the next- 
generation detectors such as GENIUS and XENON. To 
compute the spin-independent nucleon-neutralino elas¬ 
tic cross section the low-energy running values of the 
soft supersymmetry-breaking terms at the scale M^, as 
obtained b y S uSpect, were passed to the Fortran code 
DarkSUSY |^. The result of this analysis is the spin- 
independent elastic cross section cr^p^. 

The ability of a given detector to observe a signal from 
relic neutralino elastic scattering is a function of both the 
interaction cross section cr^p^ and the number density of 
LSPs in the local halo. Rather than plot event rates per 
kilogram of target per year, we have chosen to follow the 
custom in theoretical work on direct detection and 
plot an effective cross section as a function of the neu¬ 
tralino mass rn^o. This involves a rescaling of the cross 
section, normalized to the case where the relic density in 
the cosmos is the cosmologically preferred value. In our 
case we have chosen the procedure outlined by Bottino et 
al. which is to normalize the effective cross sections 
whenever the relic density, as computed by micrOMEGAs, 
is below (U^h^)caic = 0.1 by the ratio ^ = (U^h^)caic/0.1: 


^ These regions may also be accessjWe to the proposed liquid- 
Xenon based ZEPLIN 4 detector |53[ whose reach is somewhat 
less than the GENIUS reach shown here. 


eff 

XP 


t ^calc 




calc \ _calc 


( 6 ) 


This procedure allows us to consider cases where the relic 
density of neutralinos is far below the amount needed to 
provide all of the cold dark matter of the cosmos. We feel 
such an approach is important because (a) it represents 
a logical possibility given that other sources of cold dark 
matter can be contemplated, (b) is a more likely outcome 
in cases where r < 1 than in the case of mSUGRA, and 
(c) points in parameter space incapable of accounting 
for all the needed cold dark matter are often points for 
which direct detection rates are relatively high. Having 
performed this rescaling we can now display model points 
together with the expected reach in cross section for var¬ 
ious dark matter detectors. This is the manner in which 
the dashed detection reach contours were computed for 
the figures in Section |. 

A more common method for displaying the parame¬ 
ter space reach of a particular detector is the manner 
of Figure In this figure we have sampled the pa¬ 
rameter space of Figures and that of Figures in 
intervals of 25 GeV for toq and Mi. We then eliminate 
those points that failed to satisfy the EWSB and neu¬ 
tral LSP constraints, the LEP bounds on the chargino 
and lightest CP-even Higgs masses, or had too large a 
relic density U^h^ > 0.3. The calculated elastic cross 
section is then rescaled, if necessary, to produce the ef¬ 
fective cross section cr®® which is then plotted versus 
physical LSP mass. This procedure was performed for 
tan/3 = 5 (the top panels), tan/3 = 35 (middle panels) 
and tan /3 = 50 (lower panels). The left panels labeled (a) 
enforce M 3 = Mi while the right panels labeled (b) en¬ 
force M 3 = M 2 . The dark points are those of mSUGRA 
in which r = 1.0 while the lighter points include both 
the case r = 0.60 and r = 0.65 to give some idea of the 
effect of the gaugino mass nonuniversality on the typical 
effective nucleon-neutralino elastic cross sections. 

We should note that the reach of current and near¬ 
future direct detection experiments does not extend 
to the effective cross section range displayed in these 
plots. This is in accordance with the usual result for 
models with unified scalar masses at some large initial 
scale, such as mSUGRA, where the ^ term tends to be 
large com par ed to gaugino masses at the electroweak 
scale 1^, Q 1^. Thus the LSP in all of the cases we 
consider has a relatively small Higgsino content. Also fa¬ 
miliar from the unified mSUGRA case is the observation 
that the typical size of the elastic scattering cross section 
increases with higher tan /3. When the universality rela¬ 
tion between M 2 and Mi is relaxed (while maintaining 
M 3 = Ml), as in the left panels of Figure 10, the in¬ 
creased coannihilation cross section generally depresses 
the effective cross section somewhat. 

In the right panels of Figure ^ we tie the gluino mass 
to M 2 . In this case the lighter gluino results in lighter 
squarks at the weak scale, as well as a lighter Higgs sec- 
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FIG. 10: Effective elastic cross sections and detector 
reach. The rescaled, effective spin-independent cross section 
(7p^ is plotted versus the physical LSP mass m.^o for tan (3 = 5 
(top panels), tan/3 = 35 (middle panels) and tan/3 = 50 (bot¬ 
tom panels). The left pan els labeled (a) enforce the relation 
Ms = Ml of Section IA while the right panels labeled (b) 
enforce the relation Ms = M 2 of Section IB. Dark (black) 
points are those of mSUGRA where r = 1. The lighter (green) 
points are those in the region 0.60 < r < 0.65. The estimation 
of the GENIUS and XENON detector reaches was obtained 
via the dark matter plotting program of Ref. M. 


tor though the smaller value of /i that results. This dra¬ 
matically increases the importance of those squark and 
Higgs exchange diagrams that contribute to the spin- 
independent elastic cross section for neutralinos on nu¬ 
cleons. The cross section in this paradigm is augmented 
by one to two orders of magnitude: in particular we note 
that at large tan/3 nearly the entire range of parameter 
space shown will be covered by the proposed XENON ex¬ 
periment. This marked increase in the elastic cross sec¬ 
tion is consistent with a wide array of previous studies. 
The reason is clear: any physics which reduces the size 
of the low-energy squark and Higgs masses - whether 


directly through modifying the soft terms or indirectly 
through RG effects - will increase the importance of those 
squark and Higgs exchange diagrams that contribute to 
the spin-independent elastic cross section for neutralinos 
on nucleons. This includes the gaugino mass nonuniver¬ 
salities that lower M 3 ||], beginning the RG evolution 
from a lower energy scale [pOl, or studying models with 
non-universal scalar masses ^ . 


III. FINE-TUNING AND SENSITIVITY 
ANALYSIS 

The question of fine-tuning in a given parameter, or set 
of parameters, comes down to asking how natural it is to 
end up with the parameter at a certain value. Such a de¬ 
termination should ultimately be done in the context of 
some underlying theory, such as string theory. Without 
the luxury of such a fundamental theory we are forced to 
rely on constructed notions of fine-tuning. The area of 
fine-tuning has attracted many different definitions and 
interpretations.® We are not especially fond of any par¬ 
ticular measure of tuning. Nonetheless, we think that a 
firm, standard definition should be used when compar¬ 
ing models and regions of parameter space. There are 
two common measures of fine-tuning employed when dis¬ 
cussing relic densities. 

First is a frequent, though soft, definition of fine-tuning 
we might call the “optical” or “likelihood” method. Here 
one simply looks at a plot of parameter space and “eye¬ 
balls” whether or not the allowed parameter space is fine- 
tuned by estimating the relative sizes of the preferred 
volume of parameter space to the entire volume other¬ 
wise allowed in the theo^. For instance, looking at the 
tan /3 = 5 plot in Figure Q for mSUGRA (the r = 1 limit 
for the rSUGRA models), one notices that the relic den¬ 
sity depends on a very specific ratio of toq to Mi = Mi. 
Since only a very thin sliver of parameter space meets 
the relic density constraint, one would conclude, using 
the likelihood method, that dark matter is fine-tuned in 
mSUGRA for tan/3 = 5. As we will see below, this is 
not necessarily the case when more robust measures of 
fine-tuning are employed. 

The optical method can be used on the rSUGRA pa¬ 
rameter space to come to any desired conclusion about 
fine-tuning. For instance, looking at the Mi vs. mp plot 
for r = 0.6 in Figure |^, the preferred dark matter region 
does not look especially tuned: a wide region exists that 
meets the relic density constraint. However, the opposite 
conclusion is reached if one looks instead at the Mi vs. 
r plot for tan /3 = 5 in Figure |^. As noted earlier, on 
this plot the parameter space looks quite narrow and is 
obviously rather dependent on r. 

To eliminate this ambiguity, in this paper we choose 


See, for example, the discussion of this point in Ref. 
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the second and slightly less subjective definition com¬ 
monly employed, which is based on a “sensitivity” mea¬ 
sure 1 ^, 1 ^ as applied to neutralino relic densities by 
Ellis and Olive ||] 


_ d (O^h^) 
~ Oxh 2 da 


The ai in represent parameters that are to be varied 
in the theory, typically soft breaking parameters though 
sometimes also the top and bottom Yukawa couplings 
or masses. Even though this definition is not without 
limitations [^, we find it at least provides a common 
language which makes comparisons easier. The overall 
level of fine-tuning then is given by 


TABLE I: Relic Density Sensitivity for Preferred Points in 
mSUGRA Parameter Space. 


Ml 

mo 

tan/3 

[> 

1100 

235 

5 

13.0 

1500 

343 

5 

15.7 

1100 

390 

35 

31.8 

1100 

550 

50 

78.3 

1100 

950 

50 

124.1 

1600 

740 

50 

44.2 

1600 

1500 

50 

131.9 


TABLE II: Relic Density Sensitivity for Preferred Points in 
rSUGRA Parameter Space. 




( 8 ) 


We will compute two different Aq ’s for certain sample 
points in the generalized rSUGRA parameter space. In 
the first we compute An for the case where the ai belong 
to the set {r, mo, Mi, tan/3, rrit, rrib}. Since we will 
restrict ourselves to one of the gluino mass relations of 
Sections lA and IB, this parameter set represents the 


maximal sensitivity measure for a given point. But we 
will also consider the case where a more fundamental 
theory predicts one specific value of r, so we will also 
compute the quantity Aq , where the ai now belong to the 
smaller set {mo. Mi, tan/3, mt, mb}. For example, the 
variation of r in mSUGRA has no meaning, since r = 1 
is one of the defining characteristics of this model.® Thus 
we will only compute the quantity A^ for the mSUGRA 
scenario as a baseline for comparison with the expanded 
parameter space of rSUGRA. 

We begin our sensitivity analysis by considering the 
mSUGRA case as a baseline. In Table | we provide 
the sensitivity parameter A^ for a sample of points 
which yield the cosmologically preferred relic density (the 
shaded areas in the left panels of Figures II and 1 ^). 
These values are in very good agreement with those cal¬ 
culated in Ref. ^for similar points in the parameter space. 
The sensitivity is greatest for large tan/3, which is ex¬ 
pected due to the presence of the resonant annihilation 
area which is itself sensitive to tan/3 and m^°*®. 

Let us compare these typical values with those of var¬ 
ious points in the rSUGRA parameter space. Table || 
looks at the sensitivity parameters for a set of points un¬ 
der three different gaugino mass relations: M 3 = Mi, 
M 3 = M 2 , and M 2 = Ml. The first two represent the 
case of Sections IA and 1B, respectively, while the third 


® Of course the most well-studied unified model to lay claim to 
a string theory origin - namely that of the dilaton domination 
model at tree-level — predicts even more constraints among the 
various soft parameters of the mSUGRA paradigm. 



Point 


M3 

= Ml 

M3 

= M2 

M2 

= Ml 

Ml 

mo 

tan /3 

r 

An 


An 


An 

A^ 

900 

350 

5 

0.60 

50.6 

3.3 

39.2 

15.4 

6.6 

4.5 

1100 

1500 

5 

0.60 

54.0 

3.4 

50.1 

7.7 

5.4 

5.2 

3000 

1500 

5 

0.60 

48.9 

2.5 

40.7 

16.3 

250.2 

218.6 

3000 

1375 

5 

0.60 

48.8 

2.5 

49.4 

39.4 

53.6 

47.5 

3000 

1250 

5 

0.60 

48.7 

2.5 

86.1 

73.5 

56.2 

49.7 

1100 

1500 

35 

0.60 

53.5 

3.3 

47.6 

17.2 

70.7 

69.8 

1100 

1000 

35 

0.60 

50.0 

2.5 

73.6 

68.9 

70.9 

68.8 

1600 

1500 

35 

0.60 

51.5 

2.6 

80.7 

74.4 

115.5 

112.7 

2100 

1500 

35 

0.60 

47.8 

2.7 

24.0 

6.0 

58.2 

57.8 

400 

500 

35 

0.67 

70.8 

62.9 

22.5 

22.3 

136.1 

136.0 

300 

1500 

35 

0.75 

9.9 

9.9 

20.3 

18.9 

18.0 

17.4 

1100 

750 

50 

0.60 

27.8 

23.1 

26.8 

6.0 

8.4 

7.8 

1600 

1500 

50 

0.60 

39.8 

32.7 

28.6 

7.2 

19.7 

18.8 


enforces M 3 = rMi = rM 2 . This third relation was stud¬ 
ied by the authors of Ref. || and we include it here to make 
contact with their observations. For these rSUGRA sce¬ 
narios we have chosen points which have the correct relic 
density in the M 3 = M 2 case (Figure ||), but many of the 
points also have the correct relic density in the M 3 = Mi 
case so the comparison is still instructive. 

One can see that when r represents the ratio of M 2 to 
Ml, as it does in the first two columns of Table H, its 
inclusion in the set { 0 ^} has an important effect on the 
overall degree of sensitivity. This is especially striking in 
the first column in cases where the LSP is quite massive 
(large Mi) where the correct relic density is obtained 
almost exclusively through coannihilation effects. Here 
the sensitivities are everywhere comparable to those of 
mSUGRA in Table The large impact of including r in 
the set {oi} is merely verifying that this one parameter 
effectively captures the key physics involved in the relic 
density calculation. 

In the second and third column, where the gluino mass 
is allowed to vary and is generally smaller than in the first 
column, the impact of the variable r is somewhat muted. 
As mentioned previously, the impact of M 3 is on the lower 
squark masses and the lower value of /r. In the extreme 
case of the third column of Table ||, where M 2 = Mi and 
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only the gluino mass is allowed to vary (downwards) we 
see that the impact of the relative size between M 3 and 
the other gaugino masses is minor. What is far more im¬ 
portant is the absolute magnitude of M 3 , particularly in 
points associated with the region near a Higgs resonance. 
In these examples the fine tuning is comparable to that of 
mSUGRA near a Higgs resonance and is 0(100) in size. 
This indicates that when resonant annihilation processes 
are important it is the gluino mass at the high scale that 
captures the largest fraction of the physics involved. 

Thus, in the worst case scenario for M 2 -floating 
rSUGRA, the relic density sensitivity is a factor of 4 
larger than a comparable point in mSUGRA parameter 
space. In the best case scenario, rSUGRA relic density 
sensitivities are better than comparable fine-tunings in 
mSUGRA parameter space by a factor of 4. This is true 
even ignoring the possibility that r might not be included 
in the set of variation parameters. If r is excluded, the 
relic density fine-tunings of rSUGRA compare even more 
favorably with those of mSUGRA. 

So what are we to conclude from Tables | and H? On 
the one hand we agree with the authors of Ref. ^ achiev¬ 
ing the preferred cosmological relic density requires that 
the parameter r be tuned to a specific narrow range for 
any given universal scalar mass mo, B-ino mass Mi and 
tan fj. But the same could be said for any model that 
predicts a stable relic neutralino. In mSUGRA one must 
tune tan j3 and/or the mass ratio of the stau and the light¬ 
est neutralino to a specific narrow range to get the right 
answer. In cases where the gluino is separated from the 
(still unified) B-ino and W-ino masses one must tune its 
value to a narrow range in order to achieve just the right 
value of /i at the electroweak scale. Tuning of this type is 
common to phenomenology where we lack a fundamental 
theory that can explain such relations. 

On the other hand we also agree with the authors of 
Ref. these tunings that are required of the theory are 
generally speaking not /irae-tunings. That is to say, the 
sensitivity to input parameters represented by the values 
for Aq in Tables || and || are on the whole quite mild. 
In fact the relatively small tuning parameters An are all 
the more surprising when we recall that calculating the 
relic density in the presence of coannihilations involves 
the exponential of the mass difference A of (||). If a sta¬ 
ble neutralino is a significant component of the cold dark 
matter of the cosmos (and we reiterate that it certainly 
need not be), then some specific area of the MSSM pa¬ 
rameter space will inevitably be singled out. But it is 
worth bearing in mind that the tuning of parameters im¬ 
plied in Tables | and is generally small in comparison 
with the tuning required in the electroweak symmetry¬ 
breaking sector to obtain the correct Z-boson mass - a 
physically measured quantity. 


Conclusion 


Rapidly advancing calculational tools have put the 
computation of many important observables relating to 
relic neutralinos at hand to an ever-widening community. 
We believe making sound judgements as to the region 
of supersymmetric model space favored by cosmological 
observations on the density of cold dark matter requires 
utilizing these tools in directions beyond that of unified 
models. In this work we have focused on the case of non- 
universal gaugino masses and, in particular, a paradigm 
in which the W-ino mass M 2 (and perhaps the gluino 
mass) are significantly smaller at the high-energy input 
scale than the B-ino mass Mi. Such a regime may be 
a likely outcome of effective supergravity theories based 
on the weakly-coupled heterotic string, though we have 
presented the implications of these nonuniversalities in a 
model-independent manner. 

Many complications and uncertainties exist that 
should give us pause in trying to make exact quantitative 
statements: broad theoretical questions such as whether 
there may be non-thermal LSP production mechanisms 
or whether phase transitions at late times may have gen¬ 
erated some degree of late-term inflation. There are also 
technical issues that can have a significant impact on 
one’s results, particularly at large tan (3 where the choice 
of physical pole masses for the quarks, the accuracy of 
the RG evolution employed, and the completeness of the 
corrections to Higgs masses and widths are all crucial. 

Nonetheless, it is clear that relaxing the unification 
condition on gaugino masses opens up new directions 
in the parameter space which are every bit as theory- 
motivated and phenomenologically appealing as those 
of mSUGRA. These new directions may improve the 
prospects for detecting relic neutralinos from our halo 
in future experiments, particularly the proposed liquid 
Xenon experiments. The necessary tunings involved may, 
in fact, be a tool in helping us to ascertain the correct un¬ 
derlying theory once a discovery is made (for example, se¬ 
lecting a narrow list of possible hidden sector condensing 
gauge groups in string-derived supergravity models [^). 
In this sense the work presented here on nonuniversal 
gaugino masses forms an attractive complement to the 
recent work on nonuniversal Higgs masses ^ . Both 
are well-defined departures from the baseline mSUGRA 
paradigm which can claim solid motivation from heterotic 
string inspired models. Gomputation of the relic density 
and elastic scattering cross section for relic neutralino 
LSPs is a wonderful laboratory for studying supersym¬ 
metric phenomenology in that it combines elements from 
the squark, slepton, gaugino and Higgs sectors in inter¬ 
related ways. We believe disentangling these diverse 
physics contributions is best performed with a toolkit 
that contains a variety of more general MSSM models 
with which to work. 
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